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configuration,  which  we  here  term  an  “embedded  electron  gas”,  is  approxi¬ 
mated  by  the  wide  parabolic  quantum  wells  now  being  grown  by  molecidar 
beam  epitaxy  in  the  GaAs/GaAlAs  system.  Compared  with  the  electron 
density  profile  of  the  standard  Lang-Kohn  jellium  edge,  the  boimdary  of  the 
embedded  electron  gas  is  sharper  as  a  result  of  the  continuation  of  the  jellium 
background  outside  the  electron  gas.  This  is  expected  to  cause  significant 
differences  between  the  surface  properties  of  the  embedded  electron  gas  aoid 

those  of  the  “regular”  jellium  surface  model  long  used  to  study  the  surface 
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properties  of  simple  metals.  Since  the  readily  excited  centre-of-mass  motions 
of  the  parabolic-well  electronic  density  are  now  known  to  be  insensitive  to 
subband  and  many-body  effects,  surface  properties  such  as  those  considered 
here  may  take  on  increased  significance  for  the  study  of  many-body  phenom¬ 
ena  in  such  systems. 
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ABSTRACT 

We  consider  a  nonneutral  system  consisting  of  a  broad  but  finite  slab  of  immobile 
uniform  positive  charge  interpenetrated  by  mobile  interacting  electrons  of  total  charge  less 
than  that  of  the  positive  background.  This  idealised  configuration,  which  we  here  term  an 
"embedded  electron  gas",  is  approximated  by  the  wide  parabolic  quantum  wells  now  being 
grown  by  molecular  beam  epitaxy  in  the  GaAs/GaAlAs  system.  Compared  with  the  electron 
density  profile  of  the  standard  Lang-Kohn  jellium  edge,  the  boundary  of  the  embedded 
electron  gas  is  sharper  as  a  result  of  the  continuation  of  the  jellium  background  outside  the 
electron  gas.  This  is  expected  to  cause  significant  differences  between  the  surface  properties 
of  the  embedded  electron  gas  and  those  of  the  "regular"  Jellium  surface  model  long  used  to 
study  the  surface  properties  of  simple  metals.  Since  the  readily  excited  centre-of-mass 
motions  of  the  parabolic-well  electronic  density  are  now  known  to  be  insensitive  to  subband 
and  many-body  effects,  surface  propenies  such  as  those  considered  here  may  take  on  increased 
significance  for  the  study  of  many-body  phenomena  in  such  systems. 


PACS  NUMBERS:  73.20  electronic  surface  structure,  73.20D  Quantum  wells 
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The  jellium  surface  (1.2) 

a  well  known  zeroth-order  model  for  prediction  of  the 
electronic  properties  of  simple  metal  surfaces.  It  consists  of  a  half-space  of  fixed  uniform 
positive  background  charge  occupying  the  region  z  <  0,  plus  a  neutralizing  cloud  of  mobile, 
interacting  electrons.  The  width  and  location  of  the  electron  density  falloff  region  (selvage) 
near  the  positive  background  edge  are  well  known  to  affect  surface  properties.  For  example  the 
dipole  moment  of  this  diffuse  charge  distribution  is  a  major  contributor  to  the  work 
function^^\  Other  measurable  properties,  such  as  the  dispersion  of  surface  plasmons^^’'*^  the 
existence  of  "multipole"  surface  plasmon  modes^^’^’^*®^  and  the  details  of  image  and  Van  der 
Waals  forces  on  external  particles^^*^^^  are  related  more  distandy  to  the  details  of  the  selvage. 

The  jellium  model  is  of  great  theoretical  interest  because  its  geometrical  simplicity 
(lack  of  discrete  ions)  allows  much  progress  to  be  made  on  the  difficult  many-electron 
problem.  On  the  other  hand,  its  predictions  are  hard  to  test  direcdy  by  experiments  on  metals 
because  it  is  plainly  only  a  zeroth-order  approximation  to  the  true  situation  in  which  discrete 
ions  are  present 

In  this  connection,  a  promising  development  is  the  recent  fabrication  by  Molecular 

Beam  Epitaxy  (MBE)  of  "wide  parabolic  quantum  wells"  in  the  GaAs/GaAlAs  semiconductor 
ni  12  131 

system'  ’  ’  The  A1  content  is  varied  from  layer  to  layer  in  the  z  direction  perpendicular 
to  the  epitaxial  layers  in  such  a  way  that  the  energy  Ee  of  the  conduction  band  edge  is  a 
quadratic  function  of  z.  This  mimics  the  parabolic  electrostatic  potential  well  produced  by  a 
fictitious  slab  of  uniform  positive  background,  whose  charge  density  eno  is  determined  by  the 
chosen  epitaxial  growth  profile.  Spatially  remote  doping  causes  electrons  to  fall  into  this 
parabolic  well,  where  they  have  a  long  mean  free  path.  Total  doping  levels  can  be  chosen  so 
that  the  electrons  do  not  fill  the  parabolic  well  to  its  edges,  but  rather  aggregate  in  a  central 
layer  of  number  density  approximately  no  and  a  thickness  L  less  than  the  parabolic  well 
thickness  W.  A  typical  experiment  can  achieve  no  »  IQi®  cm-3,  giving  a  Fermi  wavelength 

O 

X{  of  0(1000  A).  Because  of  the  small  band  mass  m*  s  0.07m  and  the  large  semiconductor 

O 

dielectric  constant  e  s  12,  the  effective  Bohr  radius  a^*  is  of  0(100  A)  and  the  Thomas-Fermi 

O  _ 

screening  length  is  of  0(1(XX)  A).  Thus  all  lengths  relevant  to  the  many-electron  problem  far 
exceed  the  semiconductor  lattice  spacing  a  s  5  A.  This  makes  the  GaAs/GaAlAs  system  a  far 
better  approximation  to  jellium  than  the  simple  metals  for  which  the  lattice  spacing  is 
comparable  to  the  other  lengths.  It  is  also  important  to  note  that  the  dimensionless 
interelectron  spacing  »  (4mi^)'^^ag‘^  is  around  3  so  that  the  GaAs/GaAlAs  parabolic 

well  is  a  good  model  of  an  electron  gas  of  "metallic"  density,  even  though  the  actual 
conduction  electron  density  is  much  lower  than  in  the  simple  metals. 
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This  good  analogy  with  metallic  jellium  has  led  to  the  hope  that  measurements  on  wide 
parabolic  quantum  wells  can  directly  test  the  jellium  predictions  of  many  body  theory.  It  is 
well  known,  of  course,  that  this  electron  gas  is  distinctly  finite,  at  most  a  few  Fermi 
wavelengths  wide,  in  the  z  direction  perpendicular  to  layer  growth.  This  puts  it  in  an 
interesting  crossover  region  between  two  and  three  dimensional  behaviour,  and  it  is  certainly 
necessary  in  general  to  allow  for  the  discreteness  ("subbband  structure")  of  the  energy  spacing 
for  electron  modon  in  the  z  direction.  It  has  been  shown,^^^^  however,  that  neither  the 
subband  structure  nor  many-body  effects  are  observable  in  the  simplest  resonance  experiments 
in  which  the  slab  is  excited  by  a  microwave  field  which  is  spatially  uniform  across  the 
thickness  of  the  sample. 

In  the  present  paper  we  point  out  that  the  wide  parabolic  quantum  well  differs 
systematically  from  the  regular  metallic  jellium  model  in  another  respect  which  is  not  directly 
attributable  to  discrete  subband  structure:  its  electronic  surface  density  profile  tends  to  be 
narrower,  for  a  given  value  of  the  background  density  no.  This  is  due  to  the  presence  of 
excess  unneutralised  positive  background  extending  outside  the  electron  density  edge,  so  that 
we  may  speak  of  an  "embedded  electron-gas  edge".  In  view  of  the  insensitivity^^^^  of 
uniform-excitation  experiments  to  many  body  effects,  one  is  led  to  consider  experiments  with 
nonuniform  excitation  which  will  tend  to  emphasise  the  effects  of  such  surface  details. 

The  surface  profile  of  the  "embedded"  electron  gas  depends  of  course  on  the  interi^^r 
density  no  just  as  in  the  "regular"  jellium  surface.  In  principle  the  surface  profile  could  also 
depend  in  a  crucial  way  on  the  thickness  L  of  the  electron  layer.  In  this  note  we  point  out  that 
this  is  not  in  fact  the  case  provided  that  W  »  L,  i.e.;  provided  that  electron  gas  remains  well 
embedded.  Viewed  as  a  function  of  the  electron  layer  thickness  L,  the  surface  profile  settles 
down  quite  quickly  to  its  L  -•  ••  limit  once  L  exceeds  a  Fermi  wavelength  or  so.  Thus  we  can 
speak  of  a  new  standard  surface  physics  problem,  that  of  the  embedded  electron-gas  edge, 
which  depends  only  on  the  background  density  no  and  otherwise  has  a  considerable  degree  of 
univenality. 

While  these  facts  could  be  gleaned  in  principle  by  examination  of  published  density 

profiles^^^*^^^  f(x  wide  GaAlAs/GaAs  parabolic  w^lls,  it  appean  that  this  new  edge  profile 

has  not  been  systematically  compared  with  that  of  the  "regular"  jellium  surface  model  intended 

(\  2'\ 

for  the  surfaces  of  the  simple  metals.  ’  ^  To  achieve  this  comparison  we  now  present  two 
complementary  sets  of  zero-temperature  self  consistent  Local  Density  Functional  calculations. 
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The  first  calculation  was  performed  simply  by  adding  an  infinitely  thick  layer  of 
positive  background  charge  outside  a  regular  semi  infinite  jellium.  That  is,  an  extra  external 
potential 


AV®^*(z)  =  0(z).2OT<,eV  (1) 

was  added  to  an  existing  code  intended  to  calculate  the  regular  jellium  surface  profile.  The 

self  consistent  Kohn  Sham  equations^  for  this  modified  problem  were  solved  by  simple 

iteration  with  a  mixing  fraction  v  «  1  at  each  iteration,  continuing  till  the  density  profile  near 

the  surface  was  selfconsistent  to  better  than  4  significant  figures.  The  selfconsistency  away 

from  the  surface  was  better  still.  The  resulting  edge  density  profiles  can  be  regarded  as  those 

of  infinitely  thick,  deeply  embedded  electron  gas  and  are  shown  as  solid  lines  in  Figures  [1] 

and  [2],  labelled  "Embedded,  L  -» «•".  These  two  Hgures  refer  to  different  bulk  densities, 

*  * 

corresponding  to  r,  =  2.07  and  r,  =  5  respectively.  For  comparison,  surface  profiles  for  a 
regular  jellium  edge  were  run  (with  equation  (1)  replaced  by  zero):  the  results  appear  as 
dashed  lines  in  Figs.  [1]  and  [2],  labelled  "Regular  jellium  edge  (L  ••)".  The  "embedded" 
jellium  edges  are  sharper  than  their  regular  jellium  counterparts,  being  only  about  70%  as  wide 
for  the  case  r,  »  2.07.  This  is  not  entirely  unexpected  since  the  extra  positive  background 
lying  symmetrically  outside  the  each  edge  of  the  electron  layer  in  the  embedded  case  creates, 
by  Gauss's  law,  an  additional  attraction  towards  the  central  (bulk)  region  for  any  electron 
straying  outside  the  edge.  What  is  not  so  obvious,  because  of  the  self-consistent  nature  of  the 
problem,  is  that  this  attraction  is  satisfied  principally  by  sharpening  the  edge  profile  rather  than 
by  more  marginal  changes  in  the  density  over  a  wide  interior  region.  As  the  above 
calculations  had  the  bulk  density  fixed  at  a  point  deep  inside  the  jellium,  this  matter  was 
clarified  further  by  a  second  set  of  calculations  described  below. 

The  second  type  of  calculation  performed  was  for  an  electron  gas  of  finite  thickness  L 
embedded  in  an  infinitely  thick  slab  of  jellium  background.  The  bare  potential  due  to  the 
positive  background  is  strictly  parabolic,  in  contrast  to  the  external  potential  due  to  a  finite 
thickness  of  positive  charge,  which  becomes  linear  outside  the  positive  slab:  see  Figure  3.  As 
well  as  the  background  density  no,  the  second  predetermined  quantity  for  each  run  was  the 
total  number  of  electrons  per  unit  area,  and  not  the  Fermi  energy  (nor  the  central  electron 
number  density,  in  contrast  to  the  first  set  of  calculations).  Within  the  constant  effective  mass 
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approximation  the  selfconsistent  Kohn-Sham  eigenfunctions  'F  for  an  area  A  of  embedded 
jellium  are  discrete  in  the  z  direction  perpendicular  to  the  surface  and  can  be  written 

=  A‘^^exp(i[kjj.r])  v^j(z).  (2) 

with  Kohn-Sham  eigenvalue 


2 


(3) 


Here 


Y^iz)  -» 0  as  z 


and 

V«ff(z)  =  «p(z)  +  /i,c(z). 
The  electrostatic  (Hartree)  potential  (p  satisfies 

^  a  4xe(n(z)-no), 


(4) 

(5) 

(6) 


and  the  exchange-correlation  potential 

> 

^*c(z)  *  |i(nexc)  I  „  =  n(z)  (7) 


was  taken  for  simplicity  from  the  Wigner  interpolation  approximation^^^  to  the 
exchange-correlation  energy  e^c  per  electron  of  a  uniform  electron  gas. 

The  electron  number  density  is  obtained  by  summing  over  ky  and  j, 

n(z)  - 

€.<  M 


(8) 
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where  the  total  number  of  electrons  per  unit  area 

«occ 

N  -  /n(2)  d2  >  (2*r(k^o„-2iii-|lr2  J  ep  (9) 

j-1 

takes  a  predetermined  value  conveniendy  defined  by  specifying  the  effective  electron  layer 
thickness  L  =  N/no  •  The  chemical  potential  (Fermi  energy) 

(10) 


and  the  number  of  occupied  ID  levels  (subbands)  are  selfconsistendy  determined  by 
increasing  until  the  value  of  fi  determined  via  (9)  and  (10)  satisfies  the  condition 


occ 


+r 


(11) 


The  determination  of  the  and  {v^j(2)}  for  a  given  effective  potential  V«ff(z)  is  a  shooting 
problem  involving  numerical  solution  of  the  ID  Schrodinger  equation  (4).  The  {^}  and  { yj) 
then  determine  a  new  density  n(z)  via  equations  (8)-(ll)  and  hence  a  new  V*®  from  equations 
(5)-(7). 

For  each  desired  value  of  background  density  no  =  3/(4;r  [rs*ao*]  )  (corresponding  to  a 
predetermined  A1  concentration  profile  of  the  epitaxially  grown  parabolic  quantum  well)  and 
each  desired  degree  of  well  filling  N  »  n^L  (corresponding  to  a  predetermined  value  of  the 
total  remote  n-*-  dopant  concentation  per  unit  area),  the  above  equations  were  iterated  until  the 
electron  density  near  the  edge  was  consistent  to  4  significant  figures.  For  stability  only  a  small 
fraction  v<l  of  the  new  density  and  potential  were  admixed  with  the  old  values  for  the  start  of 
each  new  iteration.  For  thick  election  layen  v  had  to  be  kept  small  but  this  condition  became 
less  stringent  for  naziow  layen  (small  L). 

The  complete  results  of  such  finite-layer  calculations  for  various  values  of  the  effective 
layer  thickness  L  are  shown  in  Fig  [4]  for  r,*  »  2.07  and  in  Fig.  [S]  for  r,*  «  5.00.  The 
election  density  n(z)  is  plotted,  where  the  coordinate  z  in  the  surface-perpendicular  direction 
has  been  measured  from  the  right-hand  edge  of  the  nominal  election  slab  (i.e.  from  a  point  L/2 
to  the  right  of  the  symmetry  point).  This  demonstrates  the  similarity  of  the  electron  density 
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edge  profiles  for  a  range  of  thicknesses  L.  This  is  particularly  so  for  dense  gases  (here  r*  = 
2.07)  where  the  Friedel  oscillations  are  least  pronounced.  The  same  set  of  finite-thickness  data 
are  also  presented  (lines  with  symbols)  in  figures  [1]  and  [2],  on  an  expanded  horizontal  scale 
to  emphasise  edge  details.  These  figures  show  that,  for  a  given  background  "density”  no,  any 
one  of  the  embedded-edge  profiles  is  substantially  sharper  (falls  off  over  a  shorter  distance) 
than  the  standard  jellium  edge  profile  for  the  same  background  density.  This  is  most 
pronounced  at  high  background  density  (low  r,).  For  example,  at  r,  ~  2  the  width  of  the 
embedded  edge  is  only  about  70%  of  that  of  the  standard  r,  =  2  jellium  edge^^\  Also,  at  least 
at  r,  =  2.07,  the  stability  of  the  edge  profile  is  truly  remarkable,  the  low-density  outer  region 
remaining  the  same  as  the  thick-layer  profile  even  when  the  layer  thickness  L  is  much  less 
than  a  Fermi  wavelength. 

In  summary,  we  have  shown  that 

(i)  there  exists  an  embedded  electron-gas  surface  density  profile,  dependent  on  the 
positive  background  density  no  but  otherwise  universal  for  all  well-embedded  electron  gas 
layers  except  the  very  thiimest;  and 

(ii)  this  embedded  surface  density  profile  is  narrower  than  the  corresponding  profile  of 
a  "regular"  jellium  surface. 

Although  we  have  only  investigated  the  equilibrium  electron  density  profile,  it  seems 
likely  that  other  properties  of  embedded  electron-gas  surface  will  differ  from  those  of  the 
standard  jellium  edge.  For  example,  it  has  been  shown^^^^  using  exact  operator  algebraic 
techniques  that  there  exists  an  oscillation  mode  of  an  embedded  electron-gas  slab  whose 
frequency  is  exactly  the  bulk  plasma  frequency:  this  is  not  the  case  for  a  conventional  jellium 
slab.  While  this  "centre-of-mass  motion"  mode  is  apparently^^^^  the  only  mode  visible  under 
direct  infrared  excitation,  the  use  of  gratings  of  sufficiently  short  period,  or  other  forms  of 
coupler  which  introduce  spatially  nonuniform  fields,  should  reveal  further  collective  modes, 
for  example  compressional  modes  in  which  the  centre  of  mass  is  stationary.  Unlike  the 
centre-of-mass  mode^^^^  some  of  these  modes  should  depend  on  exchange  and  correlation 
properties  and  so  lead  to  experimental  insights  on  the  jellium  problem.  It  seems  likely  that  the 
existence  and  dispersion  of  these  other  modes,  particularly  any  surface  modes^^\  will  also 
reflect  the  difference  between  the  embedded  and  the  regular  electron  gas. 
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RGURE  CAPTIONS 

FIG  1:  Details  of  the  edge  region  of  the  electron  number  density  profile  n(z)  for  a 

relatively  dense  positive  background  (r/  =  2.07).  The  heavy  dashed  curve  represents  the  result 
for  "regular”  jellium  in  which  the  positive  background  stops  abruptly  at  z  »  0.  All  other 
curves  are  for  "embedded"  electron  gases  of  various  effective  thicknesses  L,  where 
backgground  extends  to  -fintinity  in  the  z  direction.  The  significant  feature  is  the  difference 

between  the  "regular"  jellium  edge  profile  (heavy  dashed  curve)  and  the  "embedded"  profiles 
(all  other  curves). 

FIG  2:  As  for  Figure  1  except  that  the  number  density  is  lower  (r,*  *  5.00).  The 

differences  between  regular  and  embedded  edges  are  less  striking  here. 

FIG  3:  Comparison  of  bare  positive  background  potentials  for  a  finite  thickness  of 

regular  jellium  (dashed  curve)  and  for  a  finite  thickness  of  "embedded  electron  gas".  The 
region  of  positive  background  is  shown  by  the  heavy  line  on  the  horizontal  axis,  for  the  regular 
jellium  case,  whereas  the  background  extends  to  ±  ••  in  the  embedded  case.  The  regular 
jellium  background  potential  grows  only  linearly  outside  the  background  region  whereas  the 
growth  continues  parabolically,  causing  enhanced  electron  confinement  properties,  in  the 
embedded  case. 

FIG  4:  Comparison  of  electron  density  profiles  for  embedded  electron  gases  of  various 

effective  thicknesses  L  =  N/no,  for  a  dense  positive  background  (r,*  =  2.07).  The  significant 
feanire  is  that  all  these  gases  have  a  closely  similar  edge  profile,  regardless  of  thickness. 

FIG  5:  As  in  Rg.  4,  but  fbr  a  less  dense  background  (r,*  *  5.00).  The  surface  part  of 

the  profile  is  closely  similar  for  all  but  the  thinnest  election  layer. 
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